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All types of networks arise as intricate combinations of dyadic building blocks formed by pairs of 
vertices. In directed networks, the dyadic patterns are entirely determined by reciprocity, i.e. the 
tendency to form, or to avoid, mutual links. Reciprocity has dramatic effects on every networks 
dynamical processes and the emergence of structures like motifs and communities. The binary 
reciprocity has been extensively studied: that of weighted networks is still poorly understood. We 
introduce a general approach to it, by defining quantities capturing the observed patterns (from 
dyad-specific to vertex-specific and network-wide) and introducing analytically solved models (Ex- 
ponential Random Graphs- type). Counter-intuitively, the previous reciprocity measures based on 
the similarity of the mutual links-weights are uninformative. By contrast, our measures can classify 
different weighted networks, track the temporal evolution of a networks reciprocity, identify pat- 
terns. We show that in some networks the local reciprocity structure can be inferred from the global 
one. 



PACS numbers: Valid PACS appear here 
I. INTRODUCTION 



The study of link reciprocity in binary directed net- 
works [2 12] , or the tendency of vertex pairs to form mu- 
tual connections, has received an increasing attention in 
recent years [3HT2 ] IXd ] \16 \ 156] . Among other things, reci- 
procity has been shown to be crucial in order to classify 
[3] and model [3] directed networks, understand the ef- 
fects of network structure on dynamical processes (e.g. 
diffusion or percolation processes [SH2]), explain patterns 
of growth in out-of-equilibrium networks (as in the case 
of the Wikipedia or the World Trade Web QJJ), 
and study the onset of higher-order structures such as 
correlations (TTJ [12] and triadic motifs [331 HH [TBI ES] • 
In networks that aggregate temporal information such as 
c-mail or phone-call networks, reciprocity also provides 
a measure of the simplest feed-back process occurring in 
the network, i.e. the tendency of a vertex to respond to 
another vertex stimulus. Finally, reciprocity quantifies 
the information loss determined by projecting a directed 
network into an undirected one: if the reciprocity of the 
original network is maximum, the full directed informa- 
tion can be retrieved from the undirected projection; on 
the other hand, no reciprocity implies a maximum uncer- 
tainty about the directionality of the original links that 
have been converted into undirected ones [3]. In partic- 
ular intermediate cases, significant directed information 



can be retrieved from an undirected projection using the 
knowledge of reciprocity [S]. In general, reciprocity is 
the main quantity characterizing the possible dyadic pat- 
terns, i.e. the possible types of connections between two 
vertices. It is the combination of several such dyadic in- 
teractions that in the end produces the large-scale struc- 
ture of any network. Therefore, in order to adequately 
understand or engineer the indirect effects triggered by 
the combination of many local interactions, a complete 
characterization of patterns observed at the dyadic level, 
and hence of the reciprocity structure, is necessary. 

While the reciprocity of binary networks has been stud- 
ied extensively, that of weighted networks has received 
much less attention [T7H2T)] . Obviously, this is not be- 
cause reciprocity has less important effects in weighted 
networks. Actually, the opposite is true, due to the 
remarkably larger (actually infinite) number of dyadic 
structural possibilities in the weighted case, and the 
larger potential diversity and strength of their conse- 
quences. Rather, it is this richer phenomenology at the 
dyadic level that has long prevented or complicated a 
consistent interpretation of any measure of reciprocity 
of weighted networks. While in a binary graph it is 
straightforward to say that a link from vertex i to ver- 
tex j is reciprocated if the link from j to i is also 
there, in a weighted network there are clear complica- 
tions. Given a link of weight Wij > from vertex i to 
vertex j, how can we assess, in terms of the mutual link 



of weight Wji, whether the interaction is reciprocated? 
While Wji — (no link from j to i) clearly signals the 
absence of reciprocation, what about a value Wji > 
but such that Wji -C Wijl This complication has gener- 
ally led to two approaches to the study of directionality 
in weighted networks: one assuming (either explicitly or 
implicitly) that perfect reciprocity corresponds to sym- 
metric weights (wij = Wji) [TTJ HH [5D] , and one looking 
for deviations from such symmetry by studying net flows 
(or imbalances), defined as toy — Wji [21]. In the lat- 
ter approach, significant information about the original 
weights, including their reciprocity, is lost: the trans- 
formed links are by construction never reciprocated, even 
if the original network produces the same results as any 
other network where w[j = lOjj+Ajj and = Wji + Aij, 
which instead has a very different reciprocity structure 
(possibly also across different pairs of vertices). Since 
Aij is arbitrary, this approach is also inconsistent with 
the first one, as it cannot distinguish networks that have 
very different symmetry properties. In particular, max- 
imally asymmetric (i.e. Ay = —Wji, implying w'j, L = 
whenever w[j > 0) and maximally symmetric networks 
(i.e. A VJ > Wij + Wji, implying w'^ 10^), which are 
treated as opposite in the first approach, are indistin- 
guishable in the second one. 

In addition to the above limitations, it has become in- 
creasingly clear that the heterogeneity of vertices, which 
in weighted networks is primarily reflected into a gener- 
ally very broad distribution of the strength (total weight 
of the links entering or exiting a vertex |22j). must be 
taken into account in order to build an adequate null 
model of a network [2"3Tf3"4"] . Indeed, the different intrin- 
sic tendencies of individual vertices to establish and/or 
strengthen connections have a strong impact on many 
other structural properties, and the reciprocity is no ex- 
ception. It is therefore important to account for such ir- 
reducible heterogeneity by treating local properties such 
as the strength (or the degree in the binary case) as 
constraints defining a null model for the network. This 
means that any real network should be compared with 
a proper benchmark defined as an ensemble of random 
networks which match the constraints that have been en- 
forced. In principle, this further complicates the anal- 
ysis of reciprocity, since such constraints are typically 
enforced through time-consuming randomization algo- 
rithms that rewire the network iteratively 1231 124) and are 
generally biased [35]. Moreover, it is not clear whether 
and how those algorithms generalize to the weighted case 
[521 133] . However, recently a fast and analytical method 
providing exact expressions characterizing both binary 
and weighted networks with structural constraints has 
been proposed [56] . This allows us, for the first time, 
to have mathematical expressions characterizing the be- 
haviour of topological properties under the null model 
considered. In this paper, we extend those results in 
order to propose new mathematical definitions of reci- 
procity in the weighted case, and to evaluate their be- 
haviour exactly under various null models that introduce 
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FIG. 1. Basic decomposition of any two dyadic fluxes (in the 
example shown, Wij = 2 and Wji = 7) into a fully reciprocated 
component (iuy = 2) and a fully non-reciprocated component 
(iVij — 5, which implies = 0). 



different constraints. This allows us to develop a consis- 
tent theoretical framework for the analysis of reciprocity 
in weighted networks, and to disentangle the effects of 
the empirical heterogeneity of vertices, as well as those of 
certain symmetries of the network, from the pure effects 
of reciprocity. In addition, we are able to assess whether 
an observed asymmetry between reciprocal links is con- 
sistent with fluctuations around a balanced but noisy av- 
erage, or whether it a statistically robust signature of 
imbalance. Finally, we introduce models that success- 
fully reproduce the observed patterns by introducing ei- 
ther a correct global reciprocity level or more stringent 
constraints on the local reciprocity structure. 



II. DEFINITIONS 

We aim at introducing measures of reciprocity which 
meet three criteria simultaneously: if applied to a bi- 
nary network, they must reduce to their well-known 
unweighted counterparts; they must allow a consistent 
analysis across all structural levels, from dyad-specific 
through vertex-specific to network-wide; they must have 
a mathematically controlled behaviour under null models 
with different constraints, thus disentangling reciprocity 
from other sources of (a)symmetry. As we show in SI, 
the above criteria naturally lead to a series of quantities 
that we now define. 



A. Dyad-specific measures 

We consider a directed weighted network specified by 
the weight matrix W , where the entry w^ indicates the 
weight of the directed link from vertex i to vertex j, in- 
cluding the case w^ = indicating the absence of such 
link. For simplicity, we assume no self-loops (i.e. wu = 
Vi), as the latter carry no information about reciprocity 
(in any case, allowing for self-loops is straightforward in 
our approach). As Fig. 1 shows, we can always decom- 
pose each pair (wij,Wji) of reciprocal links into a bidi- 
rectional (fully reciprocated) interaction, plus a unidirec- 
tional (non reciprocated) interaction. 

Formally, we can define the reciprocated weight be- 
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tween i and j (the symmetric part) as 



w£ = mm[wij,Wji] = 



(1) 



and the non-reciprocated weight from i to j (the asym- 
metric part) as 



(2) 



Note that if wjt > then w~! 



0, which makes the 
unidirectionality manifest. We can also define 



(3) 



as the non-reciprocated weight from j to z, and re- 
state the unidirectionality property in terms of the fact 
that w~j and wfj cannot be both nonzero. Thus 
any dyad (Wij,Wji) can be equivalently decomposed as 
(w^j ,wjj ,Wjj). If the network is binary, all the above 
variables are either or 1 and our decomposition coin- 
cides with a well studied dyadic decomposition [3HB] ■ 



B. Vertex-specific measures 

From the above fundamental dyadic quantities it is 
possible to define reciprocity measures at the more ag- 
gregate level of vertices. We recall that the out- and 
in-strength of a vertex i are defined as the sum of the 
weights of the out-going and in-coming links respectively: 



W j 



(4) 



In analogy with the so-called degree sequence in binary 
networks, we denote the vector of values {s° ut } as the 
out-strength sequence, and the vector of values {s™} as 
the in-strength sequence. Using eqs.([lj3|, we can split 
the above quantities into their reciprocated and non- 
reciprocated contributions, as has been proposed for ver- 
tex degrees in binary networks [31 0] . We first define the 
reciprocated strength 



E 



(/•', 



which measures the overlap between the in-strength and 
the out-strength of vertex i, i.e. the portion of strength of 
that vertex which is fully reciprocated by its neighbours. 
Then we define the non-reciprocated out-strength as 



sT l - si 



and the non-reciprocated in-strength as 



»T = E w ^ = s >™ - s * 



(6) 



(7) 



The last two quantities represent the non-reciprocated 
components of s° ut and s\ n respectively, i.e. the out- 
going and in-coming fluxes which exceed the inverse 
fluxes contributed by the neighbours of vertex i. 



C. Network- wide measures 

Finally, we introduce weighted measures of reciprocity 
at the global, network-wide level. Recall that the total 
weight of the network is 

^EE^ = E s r t = E s r (8) 

i j^i i i 

Similarly, we denote the total reciprocated weight as 

Extending a common definition widely used for binary 
graphs [IH3]; we can then define the weighted reciprocity 
of a weighted network as 



W 



(10) 



If all fluxes are perfectly reciprocated (i.e. W 4 ^ = W) 
then r = 1, whereas in absence of reciprocation (i.e. 

= 0) then r — 0. In SI we discuss the difference be- 
tween our definitions and other attempts to characterize 
the reciprocity of weighted networks [T7l420| . Here it is 
enough to note that all the available approaches are based 
on the assumption that the maximum level of reciprocity 
corresponds to a symmetric network where Wij = Wji, so 
that deviations from this symmetric situation are inter- 
preted as signatures of incomplete reciprocity. In what 
follows, we show that this is actually incorrect, and that 
different strategies are required in order to characterize 
and interpret network reciprocity in the weighted case. 



III. NULL MODELS 

We start by noting that, apart from the two afore- 
mentioned extreme cases, eq.(10) cannot by itself reveal 



whether in an observed network there is a tendency to- 
wards reciprocation or anti-reciprocation. Just like its 



(5) binary counterpart, eq.(10) is informative only after a 



comparison with a null model (NM) is made, i.e. with a 
value (r)jvM expected for a network having some prop- 
erty in common (e.g. the number of vertices N and/or 
the total weight W , see SI) with the observed one. There- 
fore, (r)jvM will be a function of those properties of the 
real network. As a consequence, networks with different 
empirical values of such quantities cannot be consistently 
ranked in terms of the measured value of r, as the lat- 
ter has different reference values for different networks. 
An analogous problem is encountered in the binary case 
[3], and has been solved by introducing a transformed 
quantity that we generalize to the present setting as 



Pnm 



_ r - (r) NM 
1 - {r} nm 



(11) 



The sign of pnm is directly informative of an increased, 
with respect to the null model, tendency to reciprocate 
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(pnm > 0) or to avoid reciprocation (pnm < 0). If Pnm 
is consistent with zero (within a statistical error that we 
quantify in SI), then the observed level of reciprocity is 
compatible with what merely expected by chance under 
the null model. The introduction of pnm has two im- 
portant consequences. Firstly, networks with different 
parameters can be ranked from the most to the least re- 
ciprocal using the measured value of pnm- Secondly, and 
consequently, the reciprocity of a network that evolves in 
time can be tracked dynamically using pnm even if other 
topological properties of the network change (as is typi- 
cally the case). 

Clearly, the above considerations apply not only to the 
global quantity r, but also to the edge- and vertex-specific 
definitions we have introduced in eqs.|I][3| and ([5][7|. For 
this reason, in SI we introduce and study three important 
null models in great detail. We briefly describe these 
models below. 



A. The Weighted Random Graph model 

To start with, we consider a network model with the 
same total weight W as the real network but with no 
tendency towards or against reciprocation, i.e. a directed 
version of the Weighted Random Graph (WRG) model 
[55] . This allows us to quantify for the first time the 
baseline level of reciprocity (r)wRG expected by chance 
in a directed network with given total weight. However, 
this null model is severely limited by the fact that it 
is completely homogeneous in two respects (see SI): it 
generates networks where each vertex i has the same ex- 
pected in- and out-strength ({s\ n )wRG — {s° ut )wR.G = 
{sijwRG Vi), and moreover this value is common to all 
vertices ((si) WRG = (s) WRG Vi). 

B. The Weighted Configuration model 

A popular and more appropriate null model that pre- 
serves the observed intrinsic heterogeneity of vertices is 
one where all vertices have the same in-strength and out- 
strength as in the real network, i.e. the directed Weighted 
Configuration Model (WCM) [SU]. In such model, since 
(sT)wcm = sf 1 and (s° ut ) WC M = s° ut Vz, the two 
sources of homogeneity characterizing the WRG are both 
absent: each vertex has different values of the in-strength 
and out-strength, and these values are also heteroge- 
neously distributed across vertices. In other words, this 
model preserves the in- and out-strength sequences sep- 
arately. 

C. The Balanced Configuration model 

Another important null model that we introduce 
here for the first time is one that allows us to conclude 
whether the observed asymmetry of fluxes is consistent 



with a fluctuation around a balanced network (i.e. 
one where the net flow at each vertex is zero). This 
model, that we denote as the Balanced Configuration 
Model (BCM), is somewhat intermediate between the 
above two models, as it assumes (like the WRG) that 
the expected in- and out-strength of each vertex are 
the same, i.e. that the two observed values sf 1 and 
g out are fluctuations around a common expected value 
(si) bcm — {sf +s° ut )/2, but at the same time preserves 
(as the WCM) the strong heterogeneity of vertices (i.e. 
in general (s^bcm (sj)bcm if i ^ j). This model 
preserves the total strength s* ot = sf 1 + s° ut of each 
vertex, but not the in- and out-strength separately. 

Note that all the above null models preserve the total 
weight of the original network, i.e. (W)mm — W. How- 
ever, they do not automatically preserve the reciprocity 
(neither locally nor globally). Our aim is to understand 
whether the observed reciprocity can be simply repro- 
duced by one of the null models (and is therefore triv- 
ial), or whether it deviates systematically from the null 
expectations. In the next section we show that the latter 
is true, and that the reciprocity structure is a robust and 
novel pattern characterizing weighted networks. 



D. A unifying formalism 

As we show in SI, it is possible to characterize all 
the above null models analytically, and thus to calcu- 
late the required expected values exactly. Even if the 
final expressions are rather simple, their derivation is in 
some cases quite involved and requires further develop- 
ments of mathematical results that have appeared rela- 
tively recently in the literature [251 [Ml [55] , Moreover, 
the crucial step that fixes the values of the parameters 
of all models requires the application of a maximum- 
likelihood method that has been proposed by two of us 
only recently |56j . It is for the above reasons, we be- 
lieve, that the reciprocity of weighted networks has not 
been studied as intensively as its binary counterpart so 
far. By putting all the pieces together, we are finally 
able to approach the problem in a consistent and rigor- 
ous way. Importantly, the framework wherein our null 
models are introduced (maximum-entropy ensembles of 
weighted networks with given properties) extends to the 
weighted case, and at the same time formally unifies, 
recent randomization approaches proposed by physicists 
and well-established models of social networks introduced 
by statisticians, i.e. the so-called Exponential Random 
Graphs or p* models (see SI). While a variety of spec- 
ifications for the latter exist in the binary graph case 
[H [37J |3H] , very few results for weighted graphs are avail- 
able [3J5]. Our contribution opens the way for the intro- 
duction of more general families of Exponential Random 
Graphs for weighted networks. Indeed, besides the null 
models discussed above, we will also introduce the first 
models that correctly reproduce the observed reciprocity 
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structure, either at the global (but not necessarily local) 
level, or at the local (and consequently also global) level. 
It is worth mentioning that our approach makes use of 
exact analytical expressions, and allows to find the cor- 
rect values of the parameters both in the null models and 
in the models with reciprocity. By contrast, the common 
methods available in social network analysis to estimate 
binary Exponential Random Graphs rely on approximate 
techniques such as Markov Chain Monte Carlo or pseudo- 
likelihood approaches [2 [37] [38] . Another advantage is 
that the method we employ allows us to obtain the ex- 
pected value of any topological property mathematically, 
and in a time as short as that required in order to measure 
the same property on the original network [56] . Unlike 
other randomization approaches [531 1M] > we do not need 
to computationally generate several randomized variants 
of the original network and take (approximate, and gen- 
erally biased |35j ) sample averages over them. 



interactions. We conclude that previous measures based 
on symmetric expectations are inadequate indicators of 
reciprocity Similarly, studies of asymmetry focusing on 
the differences lUy — Wji are severely limited by the fact 
that the observed imbalances might actually be fluctua- 
tions around a zero average (wtj — Wji) = 0, irrespective 
of the level of reciprocity. Thus reciprocity and symmetry 
are two completely different structural aspects (see also 
|10j). Even more generally, in SI we show that measures 
based on the correlation or similarity between toy and 
Wji are flawed. By contrast, as we show below, our def- 
initions succeed in capturing the increased or decreased 
tendency of a network to reciprocate link weights, besides 
the baseline similarity implied (when present) by other 
conservation laws. 



B. Reciprocity rankings classify weighted networks 



IV. RESULTS AND DISCUSSION 

Comparing real data with the above null models, and 
the null models among themselves, allows us to separate 
different sources of heterogeneity observed in networks. 
This is a key step towards understanding the origin of 
the reciprocity structure of real weighted networks. 



A. Reciprocity versus symmetry 

Rather surprisingly, a first general conclusion we can 
draw is that, independently of other properties of the ob- 
served network, the symmetry of weights (i.e. w,j = Wu) 
is completely uninformative about the reciprocity struc- 
ture, for two reasons. First, in networks with broadly 
distributed strengths (as in most real-world cases) the 
attainable level of symmetry strongly depends on the in- 
and out-strengths of the end-point vertices. It is easy to 
see this by noticing that, unless s\ n — s° ut for all ver- 
tices, it becomes more and more difficult, as the hetero- 
geneity of strengths across vertices increases, to match 
all the constraints required to ensure that Wij — Wji for 
all pairs. Therefore, even networks that maximize the 
level of reciprocity, given the values of the strengths of 
all vertices, are in general not symmetric. On the other 
hand, in networks with balance of flows at the vertex 



level (s\ 



for all vertices) an average symmetry of 



weights is automatically achieved by pure chance, even 
without introducing a tendency to reciprocate. This is 
evident from our study of the WRG and BCM, which 



both yield (wi 



{w. 



for all pairs of vertices (see SI). 



In many real networks (including examples we study be- 
low), the balance of flows at the vertex level is actually 
realized, either exactly or approximately, as the result of 
conservation laws (e.g. mass or current balance) or other 
constraints. In those networks, the symmetry of weights 
should not be interpreted as a preference for reciprocated 



We now carry out an empirical analysis of several 
real weighted networks using our definitions introduced 
above. We start with the global quantities r and pnm 
defined in eqs.(10l and (11). In Table 1 we report the 



analysis of 70 biological, social and economic networks 

@Msn]. 

We find that all networks display a nontrivial weighted 
reciprocity structure (i.e. p ^ 0), which differs from that 
predicted by the 3 null models considered (WCM, BCM, 
WRG). This means that the imposed constraints can- 
not account for the observed reciprocity. Remarkably, 
we also find that networks of the same type systemati- 
cally display similar values of p: for a given choice of the 
null model, the resulting reciprocity ranking provides a 
consistent classification of networks. This indicates that 
the reciprocity structure is a robust structural feature of 
weighted networks (by contrast, note that the ranking 
based on r would create some mixing or overlap, even if 
weak, between different types of networks). 

Some networks (social networks [4TH46] and the World 
Trade Web [47]) always show a positive reciprocity, while 
others (foodwebs [50] ) always show a negative reciprocity, 
irrespective of the null model. However, other networks 
(interbank networks [48]) are classified as weakly but pos- 
itively reciprocal under the WCM, but as strongly nega- 
tively reciprocal under the BCM and the WRG. In one 
case (neural network [2H])) the estimated level of reci- 
procity can be slightly positive, negative, or even consis- 
tent with zero depending on the null model. So, in gen- 
eral different null models provide different estimates of 
reciprocity, and consequently rank the same networks dif- 
ferently. For instance, the 5 interbank networks are more 
reciprocal than the neural network under the WCM, 
while the ranking is inverted under the BCM and the 
WRG. The WCM is the most conservative model since 
it preserves most information from empirical data, so we 
choose to rank the networks in the Table using pwcm- 

Importantly, we find that all weighted rankings are 
quite different from the ranking based on a binary anal- 
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Social networks (3 nets) 














Most reciprocal 


0.75 ±0.01 


0.75 ±0.01 


0.75 ±0.01 


0.85 


± 


0.01 


Least reciprocal 


0.59 ±0.01 


0.58 ±0.02 


0.57 ±0.02 


0.78 


± 


0.01 


World Trade Web (53 nets) 














Most reciprocal 


0.59 ±0.03 


0.57 ±0.03 


0.57 ±0.04 


0.79 


± 


0.02 


Least reciprocal 


0.43 ± 0.03 


0.35 ±0.05 


0.36 ±0.05 


0.66 


± 


0.02 


Interbank networks (5 nets) 














Most reciprocal 


0.07 ±0.02 


-0.26 ±0.03 


-0.26 ±0.03 


0.37 


± 


0.01 


Least reciprocal 


0.02 ±0.01 


-0.40 ±0.02 


-0.40 ±0.02 


0.30 


± 


0.01 


Neural network (1 net) 














C. Elegans 


0.02 ±0.01 


-0.11 ±0.03 


-0.007 ±0.01 


0.08 


± 


0.01 


Foodwebs (8 nets) 














Most reciprocal 


-0.14 ±0.26 


-0.67 ±0.20 


-0.65 ±0.20 


0.17 


± 


0.02 


Least reciprocal 


-0.34 ±0.22 


-0.97 ±0.02 


-0.97 ±0.02 


0.01 


± 


0.02 



TABLE I. Reciprocity of 70 real weighted networks (see SI for a description of the data), measured using pnm under 3 null 
models (Weighted Configuration Model, Balanced Configuration Model, Weighted Random Graph), and comparison with r. 



ysis [5]. While the various snapshots of the World 
Trade Web are systematically found to be strongly and 
sometimes almost perfectly reciprocal in the binary case 
(0.68 < pug < 0.95 under the binary Random Graph 
model [3]), here we find them to be less reciprocal than 
social networks if the additional weighted information is 
taken into account. Also, while the neural network of C. 
elegans has a strong binary reciprocity (prg = 0.41 [3]), 
here we find it to have a very weak (under the WCM), 
consistent with zero (under the WRG), or even negative 
(under the BCM) weighted reciprocity. These important, 
and unpredictable a priori, differences between binary 
and weighted estimates confirm that the reciprocity of 
weighted networks is nontrivial and irreducible to that of 
simplified binary projections. 
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C. The role of node imbalances 

Exploiting the different pieces of information provided 
by the three null models is revealing. Recall that the 
two differences between the WCM and the WRG are 
node imbalance {(sf 1 ) and (s° ut ) are equal in the WRG, 
whereas they are fixed to the observed, and in general dif- 
ferent, values in the WCM) and node heterogeneity (the 
expected strenghts of all vertices are equal in the WRG, 
while distributed as broadly as the empirical values in the 
WCM). We can use the BCM as an intermediate model in 
order to understand whether these two differences play a 
comparable role in determining the observed differences 
between the values of Pwcm and Pwrg-, or whether one 
factor dominates the other. The BCM preserves node 
heterogeneity but assumes node balance by regarding the 
observed difference between the in- and out-strength of 
each vertex as a statistical fluctuation around a balanced 
average. As we show in Fig. 2, some real networks (such 
as foodwebs and the World Trade Web) indeed appear to 
display very small fluctuations around this type of node 
balance. 

In other networks (e.g. the neural network or social 



FIG. 2. In-strength s\ n versus out-strength s° ut in four 
weighted networks in increasing order of reciprocity p: a) the 
Everglades Marshes foodweb, b) the neural network of C. el- 
egans, c) the World Trade Web in the year 2000, and d) the 
social network of a fraternity at West Virginia College (note 
that the increase in reciprocity is not necessarily associated 
with an increase in symmetry). 



networks) , the relation between the in- and out-strengths 
of vertices is much more scattered or does not hold. In 
foodwebs, that can be understood as transportation net- 
works [5T] where edges represent stationary flows of en- 
ergy among species, the almost perfect equality is due 
to an approximate biomass or energy balance condition 
at each vertex. In the World Trade Web, where edges 
represent the amount of trade among world countries, an 
approximate balance of vertex flows is economically ex- 
pected because countries tend to minimize the difference 
between their total import and their total exports, i.e. 
they try to balance their payments |52j . 

As we show in SI, the balance of vertex flows implies 
that, even without introducing a tendency to reciprocate, 
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FIG. 3. Temporal evolution of the reciprocity of the World 
Trade Web during the 53 years from 1948 to 2000: a) ob- 
served value of r (blue) and its expected values (r)jvjw under 
the Weighted Configuration Model (red), the Balanced Con- 
figuration Model (green), and the Weighted Random Graph 
(orange); b) evolution of pnm under the same 3 null models 
as above. 

the expected mutual weights are equal: (toy) = (iVji). 
This implies a larger expected reciprocated weight (W**) 
in the BCM than in the WCM and yields Pwcm > 
Pbcm, as confirmed by Table 1. However, we find that 
Pbcm and pwrg are always very similar, while they can 
be very different from pwcm- This means that the em- 
pirical node imbalances, even when very weak, can have 
a major effect on the expected level of reciprocity. Sur- 
prisingly, we find that this effect is much stronger than 
that of the strikingly more pronounced node heterogene- 
ity. However, irrespectively of whether we control or not 
for this potential effect (either by interpreting the ob- 
served asymmetries as zero-average fluctuations as in the 
BCM, or by controlling for the realized imbalances as in 
the WCM) , we still find that real data are systematically 
not reproduced by the corresponding null model. There- 
fore our analysis can successfully filter out the effects of 
(either exact or approximate) flux balances or other sym- 
metries, and quantify the statistically significant level of 
residual reciprocity. Disentangling reciprocity from sym- 
metry can have unexpected and counter-intuitive impli- 
cations: the most reciprocal of the four networks (the 
social network, see Table 1) is one of the least symmet- 
ric ones (see Fig. 2c), whereas the least reciprocal of 
the four networks (the foodweb, see Table 1) is the most 
symmetric one (see Fig. 2a). 



D. Time evolution and fluctuations 

As we mentioned, the fact that using p we can consis- 
tently rank networks with different properties in terms of 
their reciprocity also implies that we can use p in order 
to track the evolution of the reciprocity of a weighted 
network over time. For this reason, in our dataset we 
have included 53 yearly snapshots of the World Trade 
Web, from year 1948 to 2000 [47 . Note that the number 
of vertices, total weight, and other properties of the net- 
work change in time [9] . In Fig. 3 we show the evolution 
of r, (r) and p under the three null models. 

The plots confirm that the trend of r is not an ade- 



quate indicator of the evolution of reciprocity, since the 
baseline expected value (r) (under any null model) also 
changes in time as a sort of moving target (Fig. 3a). 
By contrast, the evolution of p correctly filters this ef- 
fect out and confirms that, throughout the evolution of 
the World Trade Web, the empirical level of reciprocity 
is much higher than that (p = 0) expected by controlling 
for any of the constraints enforced (Fig. 3b). 

The differences among the null models lead to in- 
teresting insights. We first note that (r)wRG and 
(f) bcm are very similar, while they are different from 
(r)wCM, confirming our previous discussion. More in- 
terestingly, we also note that (r)wc'M fluctuates much 
more than (r) wrg and (r) bcm (which instead evolve 
very smoothly), and its fluctuations resemble those of 
the observed value r (see Fig. 3a). This is due to the 
fact that, while in all the snapshots of the network the 
data confirm the existence of small fluctuations around 
the identity s|™ = s° ut (like those in shown in Fig. 2c 
for the year 2000), these fluctuations have different en- 
tities in different years. This produces the fluctuations 
observed in the evolution of (r)wc'M, and partly also 
in the observed value r. In other words, the temporal 
fluctuations of the reciprocity r are in part due to fluc- 
tuations of the empirical strengths of verticex around a 
balanced flux condition s^ n ss s° ut , confirming the impor- 
tant role of node (im)balances. After controlling for the 
time- varying entity of node imbalances (using the WCM) , 
we indeed find that the amplitude of the fluctuations of 
Pwcm have a smaller amplitude than those of p bcm and 
Pwrg (see Fig. 3b). However, the fluctuations of r and 
( r )wCM do not cancel out completely, and their result- 
ing net effect (the trend of pwcm) is still significant, and 
actually indicates the strongest level of reciprocity across 
the three null models. 

Whatever the null model considered, the temporal 
trend of p shows that the weighted reciprocity changes 
considerably in time and undergoes rises, falls, and rel- 
atively stable phases. By contrast, a binary analysis 
of the same network [9l [10] detected an almost mono- 
tonic increase of the reciprocity, with a marked accelera- 
tion in the 90's. Surprisingly, we find that the weighted 
reciprocity has instead undergone a rapid decrease over 
the same decade. This counter-intuitive result (strong 
increase in binary reciprocity and strong decrease in 
weighted reciprocity) confirms that the information con- 
veyed by a weighted analysis of reciprocity is nontrivial, 
and can be very different from the binary picture. 



E. Local reciprocity structure 

We now focus on the reciprocity structure at the lo- 
cal level of vertices, i.e. on the reciprocated and non- 
reciprocated strength st*, s"f and s"* defined in eqs.([5^ 
[7]). As clear from eq.Q, this allows us to analyse how 
different vertices contribute to the overall value of W++ 
and hence to r. As for the latter, measuring the (non- 
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FIG. 4. Relationship between total (s\ ot ) and reciprocated 
(st*) or non-reciprocated (s*~ + s~*) strength in four weighted 
networks in increasing order of reciprocity p: a) the Ever- 
glades Marshes foodweb, b) the neural network of C. elegans, 
c) the World Trade Web in the year 2000, and d) the social 
network of a fraternity at West Virginia College (black: real 
data, green: Weighted Configuration Model, red: Balanced 
Configuration Model). 



) reciprocated strength alone does not allow us to assess 
whether the vertex-specific reciprocity structure is signif- 
icant, or whether it is merely a consequence of the local 
topological properties of vertices. To this end, we com- 
pare the observed value of s^, s~[~ and s~* with their 
expected values under the WCM and the BCM. Note 
that in this case the WRG is not a good choice, since it 
is completely homogeneous. By contrast, both the WCM 
and the BCM preserve the total strength s* ot = s° ut + s\ n 
of each vertex, and can thus filter out the effects of the 
observed heterogenity of vertices. In Fig. 4 we show the 
observed and expected values of the (non-)reciprocated 
strength versus the total strength s' ot for the same four 
networks we have already studied in Fig. 2, which are 
ordered by increasing reciprocity. 

For the anti-reciprocal networks with p < (the food- 
web and, under some null model, the neural network), 
the dominant and less fluctuating contribution to s* ot 

comes from the non-reciprocated strength, and therefore close node-balance relationship s° 



the separation between the two trends is exponentially 
larger than it appears in a log-log plot). This shows that 
each vertex contributes, roughly proportionally to its to- 
tal strength, to the overall anti-reciprocity of this net- 
work (W** < (W**)nm and hence pnm < 0, see Table 
1). By contrast, in the neural network (Fig. 4b) some 
vertices (mostly, but not uniquely those with large s' ot ) 
have a larger non-reciprocated strength than expected 
under the null models, while for other vertices (mostly 
those with small s\ ot ) the opposite is true. This shows 
that the weak (and nearly consistent with zero, see Ta- 
ble I) overall reciprocity of this network is the result of 
several opposite contributions of different vertices, that 
cancel each other almost completely. The World Trade 
Web (Fig. 4c) also shows a combination of deviations in 
both directions, even if in this case for the vast major- 
ity of vertices the observed reciprocated strength is larger 
than the expected one. This results in the overall posi- 
tive reciprocity of the network, but again in a such a way 
that the global information is not reflected equally into 
the local one. Finally, the social network (Fig. 4d) dis- 
plays a behaviour analogous, but opposite, to that of the 
foodweb: the observed reciprocated strength of each ver- 
tex systematically exceeds its expected value and gives 
a proportional contribution to the overall positive reci- 
procity. The behaviour of the neural network and the 
World Trade Web warns us about the fact that global 
measures of reciprocity in general fail in capturing the 
patterns realized at a deeper level. One might expect 
that, in networks where deviations are instead homo- 
geneous across vertices (as in the other two networks), 
the local reciprocity structure is entirely reproduced after 
adding an overall constraint on the global one. However, 
as we show in the next section, this is in general not the 
case. Addressing the local reciprocity structure of real 
weighted networks invariably appears to be a necessary 
step to be made. 



Note that the striking similarity between the predic- 
tions of the WCM and the BCM in the foodweb and in 
the World Trade Web is not surprising, since the very 

i sT in these two 



we choose to plot + versus s ° (Fig. 4a-b). Con- 
versely, for the positively reciprocal networks with p > 
(the World Trade Web and the social network) the dom- 
inant contribution comes from the reciprocated strength, 
so we consider versus s*°* (Fig.4c-d). 

We found very rich and diverse patterns. In all net- 
works, the selected quantity displays an approximately 
monotonic increase with s' ot . Qualitatively, this increas- 
ing trend is also reproduced by the two null models. How- 
ever, we systematically find large differences between the 
latter and real data. In the foodweb (Fig. 4a), the ob- 
served values of the non-reciprocated strength + 
are always larger than the expected values (note that 



networks (see Fig. 2a and 2c) implies a correspondingly 
tight agreement between the WCM and the BCM. How- 
ever, in the neural network and in the social network the 
similarity between the predictions of the two null mod- 
els is nontrivial, since node balance is strongly violated 
in these cases (see Fig. 2b and 2d). This shows that 
the (non-)reciprocated strength has a much more regular 
behaviour than s° ut and s\ n , and suggests that it might 
indeed represent a key quantity characterizing vertices, 
even more robustly that the in- and out-strengths them- 
selves. We will exploit this result below in order to define 
models where the (non-)reciprocated strength is enforced 
as a basic constraint. 
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V. MODELLING THE OBSERVED 
RECIPROCITY STRUCTURE 

Having shown that the reciprocity of real weighted net- 
works is very pronounced, and strikingly different from 
what expected by simply taking into account either lo- 
cal (strength sequence) or global (total weight) structural 
properties, we conclude our study by introducing a class 
of models aimed at correctly reproducing the observed 
patterns. Following the spirit of the previous analysis, 
we want to keep our discussion as general as possible by 
introducing a rigorous way to understand the effect of dif- 
ferent topological properties on each other, irrespective 
of the (necessarily different) underlying physical mech- 
anisms producing the specific values of these properties 
in different networks. To this end, rather than propos- 
ing untestable models of network formation, we expand 
the null models we have considered above by enforcing 
additional or alternative constraints on the reciprocity 
structure. Our approach therefore leads us to define the 
weighted counterparts of the binary p* models with reci- 
procity [2 [2] and their generalizations [31131 HO]- In so do- 
ing, we obtain the complete class of weighted Exponential 
Random Graphs with reciprocity. We first define three 
models that exactly reproduce, besides the observed het- 
erogeneity of the strength of vertices, the observed global 
level of reciprocity (i.e. such that W*^ = (W**) and 
W = (W), implying p — 0). Our aim is to check whether 
this is enough in order to reproduce the more detailed, 
local reciprocity structure. Then, we introduce the most 
stringent model that exactly accounts for the full reci- 
procity structure at the single vertex level, and conse- 
quently also at the global network-wide level. A detailed 
description of our models can be found in SI. 
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FIG. 5. Relationship between total (s' oi ) and reciprocated 
(st*) or non-reciprocated (sf + sj*) strength in four weighted 
networks in increasing order of reciprocity p: a) the Ever- 
glades Marshes foodweb, b) the neural network of C. elegans, 
c) the World Trade Web in the year 2000, and d) the so- 
cial network of a fraternity at West Virginia College (black: 
real data, blue: the Weighted Reciprocity Model, orange: the 
Non-reciprocated Strength Model, green: the Reciprocated 
Strength Model; all such models reproduce the global level of 
reciprocity but not necessarily the local reciprocity structure) . 



an overall p = 0). This means that, even if all vertices 
appeared to contribute evenly and proportionally to the 
global anti-reciprocity (see Fig. 4a), adding the latter as 
an overall constraint is not enough in order to capture 
the local reciprocity structure. 



A. The Weighted Reciprocity Model (WRM) 

As a first step, we consider a model where the con- 
straints are s™, s° ut (for each vertex separately, as in 
the WCM) , and additionally W** . This model is the ex- 
act analogue of the binary reciprocity model by Holland 
and Leinhardt [Tj (the prototype of all p* models). This 
model replicates the overall reciprocity r exactly. How- 
ever, as we show in SI, it has the interesting limitation 
that in general it can only reproduce networks that are 
anti-reciprocal or, more precisely, less reciprocal than the 
WCM {pwcm < 0). 

Therefore, in our analysis we can only apply it to the 
foodwebs, and in Fig. 5a we show the result on the Ever- 
glades web. For the sake of comparison with Fig. 4a, we 
plot + s~* as a function of s\ ot . We find that, quite 
surprisingly, the model does not significantly improve the 
accordance between real and expected trends produced 
by the WCM and BCM (see Fig.4a). The only difference 
with respect to the latter is that now a few vertices with 
very large s* ot lie below the expected trend, while all the 
other vertices continue to lie above it (Fig. 5a) producing 



B. The Non-reciprocated Strength Model (NSM) 

As we show in SI, a way to (potentially) generate all 
levels of reciprocity r is to consider a slightly relaxed 
model, where the constraints are s^~, s~* (for each ver- 
tex), and W** . This model is also extremely simple to 
implement since it does not actually require the param- 
eter estimation phase, and is therefore scalable to net- 
works of any size. However, it does not automatically 
reproduce the in- and out-strength sequences, therefore 
it is only appropriate for networks where sj~~ and s~* 
are the dominant contributions to s] n and s° ut respec- 
tively, so that specifying the former largely specifies the 
latter as well. This means that, even if now there are 
no mathematical restrictions, using this model is again 
only appropriate for networks with negative reciprocity 
(pwcm < 0). In Fig. 5a we show the predictions of 
this model on the foodweb. Note that, as compared to 
the previous model, now the quantity + s^* is ex- 
actly reproduced by construction (and is therefore the 
independent variable, even if it is along the ordinate), 
while s\ ot (which is now the dependent variable) is not 
reproduced, with most vertices lying above the expected 
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trend and a few dominating ones lying below it. So the 
result is even worse than before. In Fig. 5b we also 
show the performance of this model on the neural net- 
work (which actually displays pwcm ~ 0, even if it still 
has negative reciprocity under other null models, sec Ta- 
ble 1). We confirm that, even if the agreement is now 
much better, most data continue to lie either above or 
below the expected curve. These results show that the 
reciprocated strengths cannot be simply reconciled with 
the total strengths. The specification of the overall level 
of reciprocity, plus only the non-reciprocated component 
of the strength sequences, is in general not enough in 
order to fully reproduce the local observed structure. 
Note however that, as p becomes smaller and smaller, 
this model is expected to become more and more accu- 
rate, as the deviation between total strength and non- 
reciprocated strength progressively vanishes. In the ex- 
treme limit — !• 0, this model exactly reproduces all 
the strength sequences of any network, including their 
reciprocated components. 

C. The Reciprocated Strength Model (RSM) 

Next, we introduce a dual model that is appropriate 
in the opposite regime of strong positive reciprocity (i.e. 
Pwcm > 0, especially in the limit pwcm — > !-)• In this 
model, the constraints are s*^ (for each vertex) and the 
total non-reciprocated weight = W — W 4 ^ . This 

model is most appropriate for networks where is the 
dominant contribution to s\ ot . As its dual model above, 
it is very easy to implement and can be scaled to net- 
works of any size. In Fig. 5c we show the predictions of 
this model on the World Trade Web. Now sj* is the in- 
dependent variable as it is obviously always reproduced, 
while s\ ot is the dependent variable which instead turns 
out not to be reproduced for all vertices (a behaviour 
analogous, and dual, to the NSM above). In Fig. 5d 
we show the results for the social network, and in this 
case we find that the model reproduces real data remark- 
ably well. This confirms that the model is particularly 
appropriate for strongly reciprocal networks. We there- 
fore find that, as in the dual case discussed above, if 
the overall reciprocity is moderate then the constraints 
are in general not enough in order to characterize the 
local reciprocity structure. However, in networks with 
strong overall reciprocity, this model accurately (and ex- 
actly in the limit — > 0) reproduces the reciprocated 
and non-reciprocated components separately, and hence 
all the local reciprocity structure. 



D. The full model 

The above three models produce the correct level of 
global reciprocity (i.e. (r) = r or p = 0) but not nec- 
essarily the correct local reciprocity structure. In par- 
ticular, we have shown that in networks with strong (ei- 



ther positive or negative) reciprocity, the local reciprocity 
structure can be simply inferred from the global one, 
plus some information about the heterogeneity of vertices 
(some strength sequence). Conversely, in networks with 
moderate reciprocity the local patterns are irreducible 
to any overall information, and thus constitute intrinsic 
heterogeneous features. In this case, it is unavoidable 
to use a model that fully reproduces the three quanti- 
ties s^~, s~* and separately for each vertex, by treat- 
ing them as constraints. In SI we describe this model, 
that we denote as the Weighted Reciprocated Configura- 
tion Model (WRCM) in detail. Using this model, all the 
plots in Fig. 5 are automatically reproduced exactly, by 
construction. Therefore we believe that this model rep- 
resents an important starting point for future analyses of 
higher-order topological properties in weighted networks. 
In particular, we foresee two main applications. 

The first application is to the analysis of weighted 
motifs, i.e. the abundances of all topologically distinct 
subgraphs of three or four vertices [13j [14]. In the bi- 
nary case, it has been realized that such subgraphs are 
important building blocks of large networks, and that 
their abundance is not trivially explained in terms of 
the dyadic structure. This result can only be obtained 
by comparing the observed abundances with their ex- 
pectation values under a null model that separately pre- 
serves the number of reciprocated and non-reciprocated 
(in-coming and out-going) links of each vertex. In the 
weighted case, no similar analysis has been carried out 
so far, mainly because of the lack of an analogous method 
to control for the reciprocated and non-reciprocated con- 
nectivity properties separately. This is exactly what the 
WRCM we have introduced can finally allow for. 

The second application is to the problem of commu- 
nity detection [53] in weighted directed networks, i.e. 
the identification of densely connected modules of ver- 
tices. Most approaches attempt to find the partition 
of the network that maximizes the so-called modularity, 
i.e. the total difference between the observed weights of 
intra-community links and their expected values under 
the WCM. In networks where the observed reciprocity 
is not reproduced by the WCM (as all networks in the 
present study), the difference between observed and ex- 
pected weights is not necessarily due to the presence of 
community structure, as it also receives a (potentially 
strong) contribution by the reciprocity. This means that, 
in order to filter out the effects of reciprocity from com- 
munity structure, in the modularity function one should 
replace the expected values under the WCM with the 
expected values under the WRCM. 



VI. CONCLUSIONS 

Motivated by the ever-increasing gap between the 
growth of data about weighted networks and our poor 
understanding of even the simplest dyadic properties of 
the latter, we have proposed a rigorous approach to the 
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reciprocity of weighted networks. We showed that the 
available measures, being based on the similarity or cor- 
relation between the weights of reciprocal edges, fail in 
disentangling a true tendency reciprocate from the effect 
of generic symmetries such as node balance or weight 
fluctuations around a common average. We therefore in- 
troduced new definitions that consistently characterize 
the weighted reciprocity structure across different hier- 
archical scales, from the level of edges through that of 
vertices to that of the network as a whole. After control- 
ling for many factors, we found that real networks sys- 
tematically display a rich and diverse reciprocity struc- 
ture. Global rankings of reciprocity succeed in classify- 
ing different types of networks, indicating that the reci- 
procity is a robust structural property and keeps trace 
of a network's function. We also found several interest- 
ing patterns at the vertex level, and found that the local 
reciprocity structure is in most cases nontrivial, in the 
sense that it cannot be inferred from the knowledge of the 
global one, even after controlling very accurately for the 
intrinsic heterogeneity of vertices. In order to obtain our 
results, we have introduced many mathematical models 
that represent novel extensions of the class of Exponen- 
tial Random Graphs to the weighted case. Our approach 
makes use of analytical expressions and therefore, unlike 
different approximate techniques, retrieves the exact re- 
sults for the expected quantities. We believe that our 
results form an important starting point to answer many 
open questions about the effect of reciprocity on higher- 
order structural properties and on dynamical processes 
taking place on real weighted networks. 
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Appendix A: Reciprocity of binary networks 

Before considering the reciprocity of weighted net- 
works, we briefly recall the basic definitions in the binary 
case, that were originally introduced to describe the mu- 
tual relations taking place between vertex pairs [TJ |2] . 



r b = (Al) 

where L = J2t^j a ij anci L ^ = J2i^j a ij a 3 i- The above 
quantity, r b , is not independent on the link density (or 

connectance) c = N ^_^ — jv(jv-i) = ^ : on ^ ne con ~ 
trary, it can be shown that c is the expected value of 
r under the Directed Random Graph Model (DRG in 
what follows) [21 2]. In the DRG, a directed link is 
placed with probability p between any two vertices, i.e. 
{aij) dbg — Pi Vi,j (with i ^ j). This implies 

, b\ = (L") = N(N - l)p 2 = = L_ 

[ /drg ^ N(N-l)p P ~N{N-1) 

(A2) 

showing that the expected value of r b coincides with the 
fundamental parameter of this null model, and hence de- 
pends on L and N. In order to assess whether there 
is positive or negative reciprocity, one should compare 
the measured r b with its expected value (r b ) drg ■ This 
means that r b cannot be used to consistently rank net- 
works with different values of L and N, because they 
have different reference values. Also, and consequently, 
r b cannot be used to track the evolution of a network that 
changes in time, because L and/or N will also change [3]. 

2. Reciprocity as a correlation coefficient 

This is why a different definition of reciprocity was 
proposed [3] , trying to control for the time- varying prop- 
erties by means of the Pearson correlation coefficient be- 
tween the transpose elements of the adjacency matrix [? 
}■■ 

b = S^j(Qij- c )(gj»- c ) = r b -c r b - (r b ) DRG 
9 ' E^jtei-c) 2 1-c l-(r*W 

(A3) 

A symmetrical adjacency matrix (as those for binary, 
undirected networks) represents a network with the high- 
est values of r b and p (both equal to 1), whereas a fully 
asymmetrical one, with zero values mirroring unit values 
on opposite sides of the main diagonal (like a triangu- 
lar matrix), displays the lowest value, being r b = and 
p = — c/(l — c)) [3]. This meaningful definition of reci- 
procity automatically discounts density effects, i.e. the 
expectation value of r b (under the DRG). As a result, 
consistent rankings and temporal analyses become pos- 
sible in terms of p. 



1. Reciprocity as the fraction of bidirectional links 



Appendix B: Reciprocity of weighted networks 



For binary, directed networks the reciprocity is defined 
as the fraction of links having a "partner" pointing in the 
opposite direction: 



In what follows we provide additional information 
about the possible generalization of the reciprocity to 
the weighted case. 
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1. Prom binary to weighted: the first route 



By looking at eq.(A3l, it is not clear whether a gener- 



alization to the weighted case should start from the first 
term on the left (i.e. as a correlation coefficient) or from 
the last term on the right (i.e. as the normalized excess 
from a random expectation) . This ambiguity comes from 
the fact that, for weighted networks, those two terms are 
no longer equivalent (as we now show). We therefore 
start by attempting the first route, and then consider 
the second one. 

If we follow the binary recipe from left to right, we 
define the weighted reciprocity as the Pearson correlation 

T, m _^a _ Wtat ^ 



coefficient (where, as usual, w ■■ 
After some algebra, this implies 
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FIG. 6. A double-star network, with only one reciprocated 
pair of weights. 



where, in order to produce a result formally equivalent 
to eq.(A3l, we have defined the weighted analogues of r 



and c as follows: 
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(B2) 



Note that the equivalence a — c, valid for the binary 
case, no longer holds: w ^ c w . The previous expressions 
generalize the binary ones and reduce to them when sub- 
stituting the ay's in place of the Wy 's. Moreover, inter- 
estingly enough, the coefficient c w can be expressed as a 
function of the weights' distribution mean, m, and stan- 
dard deviation, s, or, in an equivalent way, as a function 
of the so-called coefficient of variation, c v — s/m, as 



rrr 



l + cl 



We could be tempted to interpret c w as the weighted 
counterpart of the binary connectance and, r as the 
weighted counterpart of eq.( |Al[ ). However, we can show 
a simple case for which the above "product-over-squares" 
definition above fails in measuring our intuitive notion of 
reciprocity. Let us consider a simple network like that in 
Fig. 1. 

If we calculate r by choosing W12 — W12 , we obtain 
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2w 12 + Ei,jY(l,2),(2,l 
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(B4) 



where the sum in the denominator includes all the 
weights different from the central ones. Now, let us imag- 
ine a second situation where u>2i = W12 + 1; the calcula- 
tions, now, would give 



2w 12 {w 12 + 1) 



+ (W12 + 1) 2 +E 



(B5) 



i,j#(l,2),(2,l) w ij 



and we would intuitively require that r 2 < ri, for ev- 
ery choice of the involved weights, because of the greater 
disparity between the two central flows. However, it can 
be shown that under certain circumstances exactly the 
opposite result is obtained, by simply changing the non- 
central weights. In fact, by choosing the latter to satisfy 
the condition 



E 



> 1012 
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(B6) 



> is ob- 



the very counter-intuitive result 

tained. This shows that eq.(Bl| is not a good choice 



(B3) 

for a weighted extension of eqlA3l. 



Before considering the alternative route, we observe 
that we could also imagine to define a slightly differ- 
ent correlation coefficient, only between the two triangu- 
lar blocks of the weighted adjacency matrix: the upper- 
diagonal one and the lower-diagonal one. This would be 
defined as 



T,i<j( W tj - W u )(Wji - Wi) 



where w u i 



E, 



) 2 Ei<jK-i-^) 



(B7) 



is the upper-diagonal mean and 



w i 



N(N- 



- N(N-l) 

is the lower-diagonal mean. Again, this 



definition has an undesirable performance. This is evi- 
dent if we imagine a matrix whose transposed entries are 
defined as and wji = Xwij (with i < j). In this case, 
we would have 
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Hi<j( w ij -w u )(Xwi 



Xw,, 



M 2 E^(A 



\w u y 



1 (B8) 



independently of the value of A! So we could arbitrarily 
rise or lower the value of A, thus making the matrix more 
and more asymmetric, without measuring this effect at 
all. Note that this circumstance is impossible in the bi- 
nary case, as all weights are forced to be either zero or 
one, and therefore the only allowed value for A is one. 

The two examples above show that correlation-based 
definitions of reciprocity, while having a satisfactory be- 
haviour in the binary case, become problematic in the 
weighted one. Unfortunately, the few attempts that have 
been proposed so far in order to characterize the reci- 
procity of weighted networks [T7H2T)] are all based on 
measures of correlation or symmetry between mutual 
weights. Later, we show that symmetry-based measures 
are also flawed. Together with our results above, this 
means that all the available measures fail in providing a 
consistent and interpretable characterizaton of the reci- 
procity of weighted networks. 



2. Prom binary to weighted: the second route 



We now consider the second route, i.e. a definition 



that starts from generalizing the last term in eq.(A3|. 
This means that we are now free to first generalize r 
in a satisfactory way, rather than as a forced effect of 
the correlation-based definition, and then calculate its 
expected value under some appropriate null model. To 
this end, we note that the binary nature of the variables 
defining r b allows us to rewrite it in a very suggestive 
way: 




FIG. 7. The quantity r as a function of A. 



Ei<j (min[w 2 



mm w; 



Wji]) 



E»<j (mm|/r, 



+ max 



(Bll) 



Another advantage of this second definition is the pos- 
sibility of mutuating from it the concept of reciprocated 
strength in the same way as the concept of reciprocated 
degree was defined: 



^ aijaji -> = ^2 mhi[wij,Wji] (B12) 



so that a very impressive definition of reciprocity can be 
given, as 



w 



(B13) 



A further feature of this quantity is its scale-invariance: 
if all the weights are multiplied by a scale factor, wtj — ¥ 



Aw, : 



does not change, as shown below: 



E 



i^ij a ij a ji 



E 



mm cu 



E 



E 



i+3 ai 3 



(B9) 



The previous relation is consistent with the intuitive 
meaning of reciprocity, as a measure of the quantity of 
mutually-exchanged flux between vertices. So we can ex- 
tend this definition to the weighted case, to obtain 



E 



E 



i+3 w n 



l 3 ' 3 Z J 
Wi 



(BIO) 



where we have defined the total reciprocated weight as 
W** = ^2i^j minjtoy, Wji]. This definition does not suf- 
fer from the same limitations of the previous one. On the 
contrary, the more the difference between mutual links, 
the less the reciprocity, because the numerator would 
not change, while the denominator would become larger. 
Note that r < 1: in fact, since we are considering pairs 
of nodes at a time, we can rewrite it as 



A E, 



(B14) 

Moreover, in the case we had a matrix with transposed 
entries defined as u>ij and Wji = XiUij (with i < j) as in 
the example considered before, we would find 



2wi 



E 4<i (A+l) 
E s<3 2w tJ 
J2 i<:j (\+l)w ZJ (A+l): 



= 1, 



2 

(A+l) 



2A 



if A > 1 
if A = 1 
if A < 1 



(B15) 



thus obtaining a continuous function with a global max- 
imum in A = 1 as it should be (see Fig. 2). 

From all this, it follows that the second route leads us 



to an appropriate weighted generalization of eq. ( A3 1 is 



Pnm 



r - (r) NM 
1 - (r)jvAf 



(B16) 
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where the expected value of r has to be computed accord- 
ing to a chosen null model (NM). Indeed, this choice 
also gives us the possibility to choose different null mod- 
els, and compare their effects on p. From r < 1, it follows 
that p < 1. 



Appendix C: Null models 

In this section we describe in detail the three null mod- 
els we considered in order to carry out our analysis. We 
adopt the formalism of Exponential Random Graphs or 
p* models, which allows to obtain maximally random en- 
sembles of networks with specified constraints. Exponen- 
tial random graphs were first introduced in social net- 
work analysis [U El [37l [38] and then recently rephrased 
within a maximum-entropy approach typical of statisti- 
cal physics ES]. We adopt the latter notation, as it is 
more practical when, rather than approaching the prob- 
lem using approximate techniques such as Markov Chain 
Monte Carlo or pseudo-likelihood [37J GH] , one can solve 
the model analytically and obtain exact results as we do 
below. 

Exponential Random Graphs are very useful when one 
needs to understand, as in our case, the expected effects 
of a given set of topological properties, C (such as the 
total weight, or the strength sequence) on the structure 
of networks. Recently, a method based on the maximum- 
likelihood principle was proposed [56] in order to fit expo- 
nential random graphs to a real-world graph G* exactly 
|56| . This method provides null models which specify 
the effects of one or more constraints on the structure of 
the particular network G*, and hence allows to empir- 
ically detect patterns in the latter, identified as devia- 
tions from the model's predictions [53]. In the method, 
maximum-entropy exponential random graphs are gen- 
erated by specifying an ensemble Q of allowed graphs, 
and by looking for the probability P(G\6) of generating 
a single graph G in the ensemble in such a way that the 
Shannon entropy 



S(6) = -yP(G\e)lnP(G\e) 



Gee 



(CI) 



is maximum, under the constraints that the probability 

is properly normalized, X^Gee = ^' ano - * na * 

the expected value 



(C) g = J2 C(G)P(G\9) 
Gee 



(C2) 



of the set C of enforced topological properties equals the 
particular value C* = C(G*) observed on the real net- 
work G*: 



In the above expressions, 9 is a vector of Langrange 
multipliers allowing to tune the value of (C)g, and 9* 
is the specific value of 9 that makes {C)g coincide with 
C* , as dictated by the maximum-likelihood principle 
[54] . The solution to the above constrained maximiza- 
tion problem is 



P{G\9* 



-H(Gj 



(C4) 



where 



H(G, 6*) = Q* ■ C{G) (C5) 
is sometimes called the graph Hamiltonian and 



Z(9* 



Gee 



-H(G,6*) 



(C6) 



is the partition function, ensuring that the probability is 
properly normalized. The above formal results translate 
into specific quantitative expectations when a particular 
choice of the constraints, C, is made. 

Once the numerical values of the Lagrange multipliers 
are found, they can be used to find the ensemble average, 
(X)* , of any topological property X of interest: 



(X)* = ]T X(G)P(G\9* 
Gee 



(C7) 



(C)$, = c* 



(C3) 



The exact computation of the expected values can be 
very diffcult. For this reason it is often necessary to rest 
on the linear approximation method even if, in what fol- 
lows, the only approximation will be that of treating the 
expected value of a ratio, as the ratio of the expected 
values: (n/d) ~ (n)/(d). 

The next subsections will be devoted to the description 
of the null modes used in the main text. 



The Directed Weighted Random Graph (WRG) 
model 



We start with the simplest case, which is the most 
direct generalization of the binary, undirected random 
graph (Erdos-Renyi) model. For an ensenble of bi- 
nary, undirected networks, it was shown [29J that, if the 
only constraint C is the total number L of links (i.e. 
H(G,9) = 6L), then the probability P(G\9) coincides 
with that of the Erdos-Renyi Random Graph Model. In 
the latter, each pair of vertices is connected with the 
same probability p, all pairs of vertices being sampled 
independently of each other. In the framework of ex- 
ponential random graphs, the probability p is simply a 
function of 9. 
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The random graph model has already been generalized 
to the undirected, weighted case [55], by considering an 
ensemble of networks with non-negative, integer-valued 
edge weights (Wij £ N, Vi,j) and imposing, as the only 
constraint, the total weight, W = J2i<j w ij- The result 
is the Undirected Weighted Random Graph model |55j . 
where each pair of vertices is still independent as in its 
binary counterpart, and connected by an edge of weight 
w with probability q(w) = p w (l — p), where p = e~ e . 

Here we introduce the directed version of the weighted 
random graph. The hamiltonian of the WRG is 



W(G*) 



W(G*) + N(N -1) 



(C13) 



To calculate p we need the expected value of r. Look- 
ing at its definition, we need the expected value of the 



minimum between and Wji'. 



(f H ) _ E^^Kj. W 3*]) 



(W) 



(W) 



(C14) 



h(g, e) = ew = eJ2' 



thus, the partition function becomes 



(C8) considering that and are independent ran- 

dom variables, the cumulative distribution for the mini- 
mum is relatively easy to calculate: 



Z{0) = J2 e- ff(G ' e) = ]T e-'Xwv 



Gee 



Gee 



+oo 



= n e - IK 1 - e ~ v (° 9 ) 

i^j Wij =0 ijij 

(provided that e~ e < 1), that is a product over the 
N(N — 1) independent random variables, identified with 
the orderd pairs of the network's N nodes. So, every 
(non-negative, integer-valued) weighted network in the 
grandcanonical ensemble has the following probability 



P(min.[wij, wji] > w) = P{wij > w)P(wji >w)= p w p w ; 

(C15) 

from this, it follows that its expected value is 



(mm[wij,Wji]} 



wrg = ^ P(mm[wij, wji] >w)= P _ 2 

w — l 

(C16) 

Now, the expected value (that is, the ensemble aver- 
age) of r, computed in correspondence of the maximum- 
likelihood parameters, can be found by using the result 
of eq. [CT3} 



P(G) 



n, 



= l[ p ^(i-p) = l[ qij ( Wij ) 

(CIO) 

by defining p = e e . Note that this parameter has 
a precise probabilistic meaning, making even more ev- 
ident the above prescription, p < 1. In fact, (fly) = 
J2w°°=o a i3li3( w io) = P = 1 - ?y(0)- According to the 
maximum-likelihood principle [54, 56J, p has to be calcu- 
lated in terms of the observed quantities, by maximizing 
the function 



In C{9) = In P(G* \9) = ^ [w*j ln( e - 9 ) + ln(l - e" 9 )] 

with respect to 9. The solution to this optimization prob- 
lem can be found by isolating 9 in the above equation 



(C12) 

and, then, by inverting eq. C12 (note that the condi- 
tion expressed by eq. C3 is satisfied because (iVij) — 



E+OO I 



II! ; 



l-pl 



(r) 



E 



WRG 



i+3 l-(p*) 2 



E 



i¥=3 1-p* 



1 +p* 



(C17) 



2. The Directed Weighted Configuration Model 
(WCM) 

This second null model is the weighted version of the 
Directed Configuration Model, fully specified by the in- 
degree and out-degree sequences [29, 56 . The weighted 
counterparts of these constraints are the in-strength and 
out-strength sequences [57] : 



H(G\9) = Y,{uiS° ut +foT) = X>i + A>ii! (C18) 

i i=/Lj 

the partition function of the WCM is 



z $) = E e ~ H{G,§) = E e E ^- (ai+ft> « = (C19) 
Gee Gee 

+oo _, 
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(where 9 = {a, /?} and provided that e~( ai+ P j ^ < 1). 
Again, it is a product over N(N — 1) independent ran- 
dom variables. The reason becomes clearer when con- 
sidering the WCM: when the contraints are local, that is 
expressable as linear combinations of the adjacency ma- 
trix elements, the partition function factorizes and the 
probability of a given configuration factorizes as well, as 
a product of the independent random variables proba- 
bility coefficients [SS]. In this case every (non- negative, 
integer-valued) weighted network in the grandcanonical 
ensemble has a probability of the following form 



P(G) 



■pa) 



H(x iyj ) w v(i - x iyj ) = n^-Ko ( c2 °) 



by defining pij = e~( ai+ @i' = e~ ai e~ l3j = XiUj. Now, 
two parameters per vertex have to be calculated in terms 
of the observed quantities: the maximum-likelihood prin- 
ciple [56] prescribes to maximize 



In £(9) = lnP(G*\9) = £ [w^ ln(x iyj ) + ln(l - x lVo )\ 

,(021) 

with respect to x and y. The solution to the optimization 
problem can be found by solving the system 



observed differences between s° ut and sf 1 are interpreted 
as statistical fluctuations around a balanced condition, 
i.e. (s° ut ) = (sf 1 ). We can start from the WCM equa- 
tions, to specify them in this particular case. The condi- 
tion s° ut ~ sf 1 implies that Xi~ yi = Zi and this reduce 
the number of equations to solve, from 2N to N: 



Si 



=> s f*( G *) = V ^ir, Vi. 



(C25) 
= ZiZj and that 
(tUji). So, under the BCM, the expected value of 



This, in turn, implies that = pji 



the minimum and of r become, respectively, 



(mm 



Wji}) 



BCM 



1-f 



(C26) 



and 



(r) 



E 



BCM 



i<j l-(Py) 2 



E 



(C27) 



which is nothing more that a simplified version of eq. 

EH 

The fundamental insight given by the BCM is that, in 
networks with node balance, i.e. where (s° ut ) = (s\ n ), 
the expected weights are symmetric: 



*r'(G*) =e j( 



(C22) 



(again, the condition expressed by eq. C3 is satisfied 
because (w, 3 ) = = i=|t = i=^t)- 

The expected value of the minimum between Wu and Wu 
can be easily found by generalizing eq. |C16| 



(min[wij, wji}) wcm 



PijPji 



1 



PijPji 



(C23) 



and the expected value of r, computed in correspondence 
of the maximum-likelihood parameters, can be found by 
using the results of eq. |C22| 



V/WCM 



ST- PijPji 

E 



(C24) 



i^j 1— Py 



3. The Balanced Configuration Model (BCM) 

In addition to the WCM, we further developed a ver- 
sion of it that is intended to model networks where the 



J 'i'}IBCM 



1 



"31/ BCM 



(C28) 



This means that, in networks where the observed differ- 
ences between s° ut and s\ n are consistent with statistical 
fluctuations around a balanced condition, one automati- 
cally expects symmetric weights, even without introduc- 
ing any tendency towards reciprocation. This shows that 
measures of reciprocity based on the symmetry of mutual 
weights necessarily receive spurious contributions from 
other sources of flow balance. This observation concludes 
our statement that previously attempted correlation- and 
symmetry-based measures [17H20] cannot properly sepa- 
rate reciprocity from other factors. 



Appendix D: From null models to true models 

The three previous null models are defined in terms of 
constraints as the total weight, the in- and out-strength 
sequences and the total-strength sequence. So, not being 
included in the list of the constraints, the reciprocity r 
and the index p were a sort of target quantities, to test 
the power of the considered null models in reproducing 
them. 

Now, we can make a step forward and include some in- 
formation about the reciprocity structure of the network. 
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1. The Weighted Reciprocity Model (WRM) 

We start by generalizing the WCM, by adding to its 
hamiltonian a sort of "global reciprocity" defined over the 
whole network, thus fixing the total number of reciprocal 
links. This means to consider, as a further constraint, 
the quantity 



(having posed x, = e~ ai , y, = e~@* and z = e~ 7 ) and, 
consequently, the probability coefficient for the generic 
configuration G is 



^(G)=n 



i<j 



(D9) 



r^minK-,^,]^^ (Dl) 

ijtj i 



Now, the maximum-likelihood principle prescribes to 
maximize 



to obtain the following Hamiltonian 



H{G\6) = + ft*i n ) + W 



(D2) 



\nP(G*\6) = Y>*H HxiVi) + < MWi) + 



+ (2t»ny Inz - In Z^ RM {6)\ (D10) 



(where 9 = {a, (3, 7}). The resolution of this null model 
is considerably semplified by considering an equivalent 
way of rewriting it, 



H(G\ 



»t* + fast + {a t + [3,+ 7 )sf ] (D3) 



with respect to x, y and z. The solution to the previous 
optimization problem can be found by solving the system 



t(G * } =E 



(a?) ff ., Vt(Dll) 



having posed, to uniform the formalism, min[w i3 , Wji] = where 
w^j and having defined 



s° ut = s- 



J 13 
Uij 



wa = + w tj 



(D4) 



Now, the most challenging calculation is about the parti- 
tion function. This can be done by rewriting the hamil- 
tonian solely in terms of the variables wjj , and w£ , 

H(G\6) = ^[(a, + P,)w^ + {a, + + (D5) 



i<j 



+ (a, + ^ + ctj + 2y)w%] 
and considering the admissible states for them: 



(D6) 



(t"y , «£> w%) = {(0, 0, N), (N+, 0, N), (0, N+, N)} 

(D7) 

where N = [0 . . . 00) and N + = [1 . . . 00). So the parti- 
tion function becomes 



Z{6) = ]T e- H ^ = Y, e- H ^ = 

Gee {wTf,wtT,wtt) 

(1 - XiXjyiyj) 



U.w 



4< J (1 - x iyj )(l - x jyi )(l - x^jy.yjz 2 ) 

^ Z WHM0 ) 

i<j 



(D8) 



x*y*0--x* y *. 



-+- 



x*x*y*y*(z*f 



12 ' e * (1 - x*y*)(l - x*x*y*y*) ' 1 - x*x*y*y*{z*y ' 

(D12) 



XjVi^-XiVj) | x*x*y*y*(z*) 2 

(1 - x*y*){l - x*x*y*y*) 1 - x*x*y*y*(z*) 2 ' 

(D13) 

(D14) 



, ^> x*x*y*y*(z*) 2 
[W ^ )S * l-x*x* jy *y*{z*r 



Now, the expected value of the minimum between Wij 
and Wji is (min^-, Wj,])Vi?M = ( w tj)e*- Even if ^ is 
possibile to write down the analytical expression of the 
expected value of r, by using it, this can be avoided, by 
considering that 



<r>5, 



WRM 



(W)g, 

this, in turn, implies that 



W*+(G*) 
W(G*) 



(D15) 



PWRM — 



I WRM 



1 - m 



WRM 



1 - r 



= 0. 



(D16) 



So, by definition, the index p is trivially reproduced by 
the WRM. 
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2. The Non-Reciprocated Strength Model (NSM) 

A second null model including the information about 
the global reciprocity structure of the network can be 
defined, starting by the WRM hamiltonian. This time, 
the imposed constraints are the in- and out-strength se- 
quences, diminished by the reciprocal strength sequence 
(see eq. D4), and the total number of reciprocal links: 



i 



Following the calculations of the WRM, the partition 
function is 



Z{8) = £ e~ H 



(G\9) 



E 



-H{G\6) 



Gee 



( w ij: 



(1 - XiXjyiyj) 



n (i-x iyj )(i-x jyi )(i-z2) 



i<3 
i<j 



(D18) 



(having posed x, = e Qi , %ji = e /3i and z = e 7 ). The 
probability coefficient for a generic configuration, G, is 



p(G) = n 



'(9) 



(D19) 



and the maximum-likelihood principle prescribes to max- 
imize 



+ {2 W ? j y\nz-\nZ% SM 0)] 



(D20) 



with respect to x, y and z. The solution to the previous 
optimization problem can be found by solving the system 



= E 

= Ei<i2«- /() 



Vi 

Vi(D21) 



where 



(z*) 2 



(l-x*y*)(l-x*x*y*y*) 



( Z *r 

(D22) 



CO 



*\2 



(l-xM)(l-x*x*y*yt) 1-W 



As for the WRM 



V/NSM 



(W* 



{W) S , 

this, in turn, implies that 



(O 2 

1 -(**)' 



W^(G*) 
W(G*) 



r; 



Pnsm 



r — (r) 



NSM 



I NSM 



l-r 



0. 



(D23) 
(D24) 

(D25) 

(D26) 



3. The Reciprocated Strength Model (RSM) 

Until now, we have defined three null models with no 
constraints about the reciprocity (the WRG, the WCM 
and the BCM) and two null models with the total number 
of reciprocal links (thus implementing a global notion of 
reciprocity), as a constraint. 

Now, we can define more refined null models, by con- 
sidering, as constraints, the local notion of reciprocity, as 



defined by eq. |B12| We start by considering the following 
hamiltonian: 



H(G\e) = aW + 6^? 



(D27) 



(where 9 = {a, 5}). The resolution of the null model de- 
scribed by this hamiltonian is, again, considerably sem- 
plified by considering the equivalent expression 

H(G\6) =Y,[aw~; +awlj + {2a + S i +5 J )w%}; (D28) 

by summing over the states defined in eq. |D7| we find the 
partition function 

Z(9) = ]T e- H ^ = Yl e~ H ^ = 
Gg9 i w ij> w tj ' w tj) 



TT y — 

11 (l-x)(l 



(1 + x) 



i<j 



(1 — x)(l — x 2 ZiZj) 



(D29) 



(having posed x = e a , Zi = e Si ) the probability coeffi- 
cient for the generic configuration G is 



p(G) = n 



x w » +w >i +2W » {z iZj ) w » (1 - x 2 z^)(l - x) 



i<3 



n 



(l + x) 

''i +w i*(ziZj) w v(l - x 2 z iZ] ){l - x) 
+ 



(D30) 
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and the likelihood function is, of course, the logarithm of 
the previous probability coefficient. The solution to this 
optimization problem prescribes to solve the following 
system 



(D31) 



where 



{x*fz*z* 



(l-{x*Y) l-{x*) 2 z*z 



"ij IS* 



( * \ 2 * * 
[X ) Zj 

1 - {x*) 2 z*z* 



(D32) 



(D33) 



This model allows to solve for the x value analitically. 
In fact, by summing eq. |D32| over the ordered pairs of 
nodes, we find 



W{G*) = 



N(N - l)x* 
(l-(^) 2 ) 



W^{G* 



(D34) 



and by solving this second-order equation w.r.t. x, and 
taking the positive solution, we have the maximum- 
likelihood estimation of this parameter. Also this model 
exactly reproduces the observed reciprocity, because 
W^(G*) = (W^)g, and W{G*) = {W) s ,. This means 

that prsm = iI(r)BSM = T=? = and tlie local quali- 
ties as the reciprocal strength sequence are now trivially 
reproduced. 



4. The Weighted Reciprocated Configuration 
Model (WRCM) 

The last step is the dcfinion of a very general null 
model, to finally include those local quantities not fixed 
by the NSM and the RSM. This implies a slight gener- 
alization of the formulas in the previous two paragraphs. 
The graph hamiltonian becomes 



H(G\9) = + P % st + Jisf) 



(D35) 



where, now, 9 = {a, f3, 7} and 



K¥*) 

with obvious meaning of the symbols (defined above). 
The partition function now becomes 



(1 - XiXjyiyj) 



(1 - Xiiij)(l - x 3 yi)(l - Xi Xj yi yj Zi Zj ) 



Y[ z whcm {9) 



7WRCM 

and the likelihood is 



(D37) 



]nP(G*\6) = M*iVi) + «)* MWi) + 

i<j 

+ (wtj)* In(ziZj) - lnZ^ flCM (0)].(D38) 

The solution to this optimization problem, with re- 
spect to x, y and z, can be found by solving the following 
system: 



3 (^)( w tj)s* 



)\ w ij)S* 



Vi 

Vi (D39) 
Vi 



where 



J ij 1 9* 



{l-x*y*){l-x*x*y*y*y 



- x*y*){l - x*x*y*y*Y 



»ij 19' 



1 



(D40) 



(D41) 



(D42) 



By the definition of the WRCM model, we not only 
recover the result that the global reciprocity is equal to 
the observed one (implying r = (r)wRCM and pwrcm 
( j. also valid for the WRM): now, all the vertex- level, 
strength sequences are exactly reproduced, impyling that 
the reciprocity is reproduced at a local level. 

The WRCM is now powerful enough to allow for the 
analysis of the weighted motifs (to understand which all 
the dyadic information has to be fixed) and for the com- 
munity detection, especially for those newtorks where the 
reciprocity plays an important role in shaping its struc- 
ture. 



Appendix E: Models: a summary 

The first six models explained in the previous sections 
can be recovered from the last and most general one (the 
WRCM) by means of simple substitutions in the graph 
hamiltonian (as shown by the following table). 
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Non-reciprocal models 


oti = ft = 7; 7i = 27, Vi 


WRCM WRG 


o-i = ft; 7i = on + ft, Vi 


WRCM -> BCM 


ctt ^ ft; 7i = Oi + ft, Vi 


WRCM -4 WCM 


Global-reciprocity models 


tt t ^ ft; 7j = on + fa + 7, Vi 


WRCM -> WRM 


a, 7^ ft; 7i = 7, Vi 


WRCM -> NSM 


Local-reciprocity models 


at = ft = a; 7i = + ft + ^i) Vi 


WRCM -> RSM 



Appendix F: The jackknife method 

The jackknife method 58, 591] is an expedient to mim- 
ick resampling and it is usually used to estimate the vari- 
ance of a given function of the population mean, f((x)) 
(being x the random variable of interest). Doing this in 
the least biased way, would imply to have a whole col- 
lection of samples. However, we observe only a single 
realization. How can compensate for the lack of such ob- 
servations? We can build a set of artificial samples by 
considering the following sets: 



S'{ = {%2, X 3 . ..X M }, 
s 2 = { x l, x 3 ■ --Xm}, 

s J M = {xi, X 2 ■ ■ .x M -i}; (Fl) 

that is a list of vectors for each of which a single observa- 
tion has been removed. Then, we calculate the so called 
jackknife averages 



$2 



J 



E 



i=£2 ■ 



SM" 



E 



i^M ■ 



M -1 M-l 
the estimates of the first two moments 



M- 1 



J\2 



Ml - 77 , A*2 = 



M 



M 



(F2) 



(F3) 



from which the estimate of the jackkinfe-standard devi- 
ation follows 



M2 



.7^2 



and, finally 



VM - 1 a 



f 



(F4) 



(F5) 



How can we implement all this for our weighted net- 
works? The quantity we are interested in is p. It is a 
function of the expected value of r, taken over the whole 
grandcanonical ensemble: (r) . By applying the jackknife 



method, we can build L artificial samples by removing 
one weight at a time. By rewriting the above formulas, 
the final estimates become 



Pnm 



r - {r) NM 



1 - r 



NAI 



Pnm 



L 

E 



(Pi, NM — PNmY 



(l-(r N u)) 2 ' 



(F6) 



(F7) 



where NM can be WRG, WCM, BCM, WRM, NSM, 
RSM, WRCM and where the sum over the index i 
means that we are summing over the realizations with 
the i-th weight removed. 



Appendix G: Description of the dataset 

In what follows a brief description of the analysed net- 
works is given. 

Interbank network. This is the network of the Ital- 
ian interbank monetary exchanges in the year 1999. 
We analysed the monthy transactions for May (N = 215, 
L = 5269), June (N = 215, L = 5229), August (N = 215, 
L = 5071), October (N = 215, L = 4712) and December 
(N = 215, L = 4685). Food webs. We analysed eight 
different food webs [50? ? ], from different ecosystems 
(lagoons, marshes, lakes, bays, estuaries, grasses, rivers), 
with a prevalence of aquatic habitats: Chesapeke Bay 
(N = 34, L = 177) and Mondego Bay (N = 46, L = 400), 
Everglades Marshes (N = 69, L = 916), Maspalomas 
Lagoon (N = 24, L = 82), Michigan Lake (N = 39, 
L = 221), St. Marks Seagrass (N = 54, L = 536), 
Crystal River Creek (N — 24, L — 125 and N = 24, 
L = 100). Neural networks. We analysed the neural 
network gS] of C. Elegans (N = 297, L = 2345). Social 
networks. We analysed three different social networks 
BK- Office, BK-Tech and BK- Fraternity. BK- 
Tech and BK-Fraternity are completely connected (that 
is, L = N(N — 1)). Bernard and Killworth (and, later, 
also with the help of Sailer), collected five sets of data 
on human interactions in bounded groups. BK- Office 
(N = 40, L = 1558) is the network of the human inter- 
actions (conversations) frequency between the employees 
of a small business-office, as recorded at time intervals 
of fifteen minutes (during two four-days periods), by an 
external observer, along a fixed route through the office 
itself. BK-Tech (N = 34, L = 1122) is the network 
of the human interactions (conversations) frequency be- 
tween collaborators in a technical research group at a 
West Virginia University, as recorded at time intervals 
of half- hour (during one five-days working week), by an 
external observer. BK-Fraternity (N — 58, L = 3306) is 
the network of the human interactions (conversations) 
frequency between the students living in a fraternity 
at a West Virginia College, as recorded by an external 
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observer at time intervals of fifteen minutes (during a 
five-days week, twenty-one hours per day) who walked 
through the public areas of the building. The World 
Trade Network. We analyse the series of yearly bi- 



lateral data on exports and imports among world coun- 
tries from the database in ref.[47], from 1948 to 2000 
(N e [82, 186] and L G [2539, 19903]). 
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